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ABSTRACT 



We present X-ray spectra of four ultraluminous infrared galaxies (ULIRGs) 
with detectable broad near-infrared emission lines produced by active galactic nu- 
clei (AGNs): Mrk 463, PKS 1345+12, IRAS 05189-2524, and IRAS 07598+6508. 
With the exception of IRAS 07598+6508, high quality X-ray spectra obtained 
with XMM or Chandra show modest (30-340 eV) equivalent widths of the 6.4 
keV Fe Ka emission line and clear signatures for absorption at a level of 4-33 x 
10 22 cm -2 for the main power-law components from the AGNs. These spectral 
properties are typical of Compton-thin AGNs, and so we estimate absorption- 
corrected 2-10 keV X-ray luminosities for the AGNs, Lx(2-10keV), using the 
Compton-thin assumption. We compare the Lx(2-10 keV) values with broad 
optical/near-infrared emission-line luminosities, and confirm a previous finding 
by Imanishi & Ueno that the Lx(2-10 keV) to broad-emission-line luminosity 
ratios in ULIRGs are systematically lower than those of moderately infrared- 
luminous type-1 AGNs. A comparison of independent energy diagnostic methods 
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suggests that the AGNs are underluminous in the 2-10 keV band with respect to 
their overall spectral energy distributions, as opposed to the broad emission lines 
being overluminous. This X-ray under-luminosity should be taken into account 
when using 2-10 keV X-ray data to investigate the energetic contribution from 
AGNs to the large infrared luminosities of ULIRGs. 

Subject headings: galaxies: active - - galaxies: nuclei - - X-rays: galaxies - 
galaxies: individual (Mrk 463, IRAS 05189-2524, PKS 1345+12, and IRAS 
07598+6508) 

1. Introduction 

The star-formation rate is one of the most important indicators of the nature of a 
galaxy; it is usually estimated from continuum or emission-line luminosities (Kennicutt 1998). 
In dusty objects, the infrared luminosity is particularly important, because the correction 
for dust extinction is generally not significant. An observed infrared luminosity can be 
converted into a star-formation rate, assuming that the luminosity is powered primarily 
by star-formation activity. Although this assumption seems to be true in the majority of 
moderately infrared-luminous galaxies with Ljr < 1O 12 L , it has been suggested that active 
galactic nuclei (AGNs) provide an important energetic contribution in ultraluminous infrared 
galaxies (ULIRGs) with L IR > 1O 12 L (Veilleux, Kim, & Sanders 1999a). Therefore, the AGN 
contribution to the observed infrared luminosity of ULIRGs needs to be properly estimated. 

To estimate this contribution, an indicator which can trace AGN power, by disentangling 
from star-formation activity, is required. Two useful indicators are currently known. The 
first is the absorption-corrected intrinsic 2-10 keV X-ray luminosity. In an AGN, strong 2-10 
keV emission is produced by the inverse-Compton process in the close vicinity of an accretion 
disk around a central supermassive black hole, whereas the 2-10 keV emission from star- 
formation is usually much weaker. The second indicator is the broad hydrogen emission lines 
in the optical and near-infrared domains, whose line widths are larger than ~1500 km s" 1 
in full- width at half maximum (FWHM). Such broad emission lines are believed to originate 
in high- velocity gas (the so-called "broad-line regions" ) around a central supermassive black 
hole in an AGN, rather than being associated with phenomena related to star formation 
(Veilleux, Sanders, & Kim 1997b). In many dust-unobscured type-1 AGNs with moderate 
infrared luminosity, both intrinsic 2-10 keV and broad emission-line luminosities have been 
derived, and a correlation between these luminosities has been found (Ward et al. 1988). If 
this correlation holds for all types of AGNs, then both of these indicators can be used to 
estimate AGN's power in a galaxy. 
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However, results using ASCA data have shown that the absorption-corrected 2-10 keV 
to broad emission- line luminosity ratios in some ULIRGs are systematically lower than those 
of dust-unobscured, moderately infrared- luminous, type-1 AGNs (Imanishi & Ueno 1999a). 
One possible explanation is the large time variability in X-ray luminosity. Some AGNs that 
were initially classified as X-ray underluminous based on X-ray data taken at one epoch 
showed normal X-ray luminosities when observed at a later epoch (Risaliti et al. 2003). A 
second possibility is the misidentification of the origin of the detected 2-10 keV emission. 
Based on modest equivalent widths of the Fe Ka line at 6.4 keV (EW 6 . 4 < several x 100 
eV) and signatures for Compton-thin absorption (N H < 10 24 cm~ 2 ) from ASCA spectra, 
absorption-corrected 2-10 keV luminosities, Lx(2-10 keV), were derived, assuming that the 
detected X-ray emission is a direct component from an AGN behind Compton-thin absorbing 
material (Awaki et al. 1991; Krolik, Madau, & Zycki 1994; Ghisellini, Haardt, & Matt 
1994). If, instead, the detected X-ray emission is a scattered/reflected component of AGN 
emission behind Compton-thick absorbing material (Nh > 10 24 cm -2 ), the intrinsic 2-10 
keV luminosity will be larger than the estimate based on the Compton-thin assumption. In 
the ULIRG, Superantennae, although no significant Fe Ka emission was seen in an ASCA 
spectrum due to limited photon statistics (Imanishi & Ueno 1999b; Pappa, Georgantopoulos, 
& Stewart 2000), a strong Fe Ka line (EW 6 . 4 ~ 1.4 keV) was reported in a higher quality 
spectrum obtained later with XMM (Braito et al. 2003a). This demonstrated that the 
detected 2-10 keV emission is produced by a scattered/reflected component, which was 
difficult to be recognized based on ASCA data. Higher quality XMM and Chandra data 
obtained at a different epoch to the ASCA era are very useful to (1) investigate whether the 
X-ray underluminosity of AGNs in ULIRGs found by Imanishi & Ueno (1999a) is caused by 
X-ray time variability and (2) better understand the origin of detected 2-10 keV emission 
using higher photon statistics at 2-10 keV and, in particular, near 6.4 keV. 

In this paper, we present XMM and Chandra X-ray spectra of ULIRGs that have de- 
tectable broad optical/near-infrared emission lines, in order to investigate the 2-10 keV to 
broad-line luminosity ratios. Throughout this paper, we adopt H = 75 km s -1 Mpc -1 . 
Luminosities derived in previous literature using different values of H were re-calculated. 
Lx(2-10 keV) denotes the absorption-corrected 2-10 keV luminosity of an AGN, unless oth- 
erwise stated. 

2. Targets 

Our targets are Mrk 463, PKS 1345+12, IRAS 05189-2524, and IRAS 07598+6508. 
Table 1 gives detailed information on these sources. 
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Mrk 463 has a double nucleus with a separation of ~4 arcsec at 1.6 yum (Surace & 
Sanders 1999) and is optically classified as a Seyfert 2 (Sanders et al. 1988). The eastern 
nucleus (Mrk 463E) is active and shows a broad (FWHM > 1500 km s _1 ) Pa/5 (A rcst = 
1.28 //m) emission line in the near-infrared (Veilleux, Goodrich, & Hill 1997a). Based on an 
ASCA spectrum, Ueno et al. (1996) estimated an EW 6 . 4 < 700 eV, N H ~ 2 x 10 23 cm" 2 , 
and found the absorption-corrected 2-10 keV luminosity to be Lx(2-10 keV) ~ 4 x 10 42 ergs 
s -1 . Mrk 463E was diagnosed to be powered by an AGN, with no detectable star-formation 
activity in the infrared (Genzel et al. 1998; Veilleux, Sanders, & Kim 1999b; Imanishi 2002). 
It was found to be underluminous in the 2-10 keV range compared with the [OIII] and broad 
optical/near-infrared line emission (Ueno et al. 1996; Imanishi & Ueno 1999a). 

PKS 1345+12 is a double-nucleus source with a separation of ~2 arcsec, with the western 
nucleus ~3 times brighter at 2 /jm than the eastern nucleus (Surace & Sanders 1999). It 
is optically classified as a Seyfert 2 (Sanders et al. 1988), and shows a broad Pact (A res t = 
1.88 jum) emission line in the near-infrared (Veilleux et al. 1997b). Imanishi & Ueno (1999b) 
estimated an EW 6 . 4 < 500 eV, N H = 1-4 x 10 22 cm" 2 , and Lx(2-10 keV) ~ 3 x 10 43 ergs s _1 , 
based on an ASCA spectrum. This source is thought to be powered by an AGN (Veilleux et 
al. 1999b). It was found to be underluminous in 2-10 keV X-rays, compared to the broad 
optical/near-infrared emission-line luminosity (Imanishi & Ueno 1999a). 

IRAS 05189—2524 is a single-nucleus source (Sanders et al. 1988), optically classified 
as a Seyfert 2 (Veilleux et al. 1999a), and shows a broad Pact emission line (Veilleux et al. 
1999b). Severgnini et al. (2001) later reported the presence of broad components of Pa/5 and 
Pact emission lines, and argued that the broad Ha (A rcst = 0.65 jum) component detected in 
the previous optical spectrum (Young et al. 1996) is a transmitted component, rather than 
a scattered/reflected component. This source was not included in Imanishi & Ueno (1999a), 
but Severgnini et al. (2001) later presented ASCA and BeppoSAX X-ray spectra which show 
EW 6 . 4 ~ 100-150 eV, N H = 4-9 x 10 22 cm" 2 , and Lx(2-10 keV) ~ 2 x 10 43 ergs s -1 . From 
the Lx(2-10 keV) value, Severgnini et al. (2001) argued that the AGN is not powerful enough 
to account for the bulk of the infrared luminosity. However, some independent diagnostics 
in the infrared suggested that this source is dominated by an AGN (Veilleux et al. 1999b; 
Soifer et al. 2000; Imanishi & Dudley 2000). 

IRAS 07598+6508 is a single-nucleus source, optically classified as a Seyfert 1 (Sanders 
et al. 1988), shows a broad Pact emission line (Taniguchi et al. 1994), and is a member of the 
class of broad absorption line (BAL) quasars (Lipari 1994; Hines & Wills 1995). No X-ray 
emission was detected by ASCA (Risaliti et al. 2000). 
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3. Observations and Data Analysis 
3.1. XMM 

Mrk 463, IRAS 05189-2524, and IRAS 07598+6508 were observed with the EPIC 
cameras (PN, MOS1, and MOS2) on board XMM, in full-frame mode. Observing details are 
shown in Table 1. Standard data analysis procedures were employed, using the XMM Science 
Analysis System (SAS version 5.4), together with standard software packages (FTOOLS 5.2 
and XSPEC 11.2). Events recorded during high-background time intervals were removed, 
and only events corresponding to pattern 0-4 for PN and pattern 0-12 for MOS were included 
in the analysis. The most recent calibration files retrieved from the XMM webpage and 
response matrices created with SAS were used for the data analysis. 

The source spectra of Mrk 463 and IRAS 05189—2524 were extracted from 30 arcsec 
radius circular regions, whereas for IRAS 07598+6508 a 22.5 arcsec circle was used. Back- 
ground spectra were extracted from nearby source-free circular regions of 60 arcsec radius. 
Spectra were binned so that each energy channel had more than 20 counts. Since the PN 
spectra have the best photon statistics, they were used for primary estimation of important 
parameters. However, to improve the statistics, final spectral fitting was performed using 
a combination of PN, MOS1, and MOS2, keeping the relative normalization free, because 
there is some uncertainty about the relative flux calibration of the three cameras. 



3.2. Chandra 

PKS 1345+12 was observed with Chandra. Table 1 summarizes the observing informa- 
tion. A 1/8 sub-array mode of the ACIS-S3 back-illuminated CCD chip was used, with a 
0.4-sec frame time. The data were reprocessed with CIAO 2.2.1 and CALDB 2.17. Only 
ASCA ("good") grades 0, 2, 3, 4, and 6 were used in the analysis. High background intervals 
were removed. 

A nuclear spectrum was extracted from a circular region with a radius of 8 arcsec. 
A background spectrum was extracted from an annular region around the nucleus, and 
subtracted from the nuclear spectrum. The count rate from the nucleus was 0.057 counts 
s _1 , and the effect of pileup was negligible. 
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4. Results 

All sources are clearly detected in the XMM or Chandra images. For PKS 1345+12, 
the high resolution Chandra image reveals that the X-ray emission comes from the western 
nucleus (Figure 1). Figure 2 shows XMM EPIC spectra of Mrk 463, IRAS 05189-2524, and 
IRAS 07598+6508, and a Chandra spectrum of PKS 1345+12. 

In the X-ray spectra of Mrk 463, PKS 1345+12, and IRAS 05189-2524, an absorbed 
hard component is clearly seen at >2 keV. This hard component is likely to originate in the 
AGN, and so is fitted with an absorbed power law. In addition to this hard component, 
there exists a soft X-ray component which dominates the emission at <2 keV. There are 
two physically plausible explanations for the origin of the soft X-ray emission: (1) a scat- 
tered/reflected component of the primary AGN power-law emission, and (2) emission from 
a starburst. In the caption of Table 2, detailed descriptions of the models used are reported 
for each source. 

For the spectral fitting of Mrk 463, a single temperature thermal component (from a 
starburst) or a second power law (a Thomson-scattered component of the primary power law 
emission from the AGN and/or X-ray binaries in a starburst) was initially included, in addi- 
tion to the main absorbed power-law component. However, significant residuals remained at 
<2 keV. We then tried model (A), a low-temperature (0.5-1 keV) thermal model + a high- 
temperature (5-10 keV) thermal component (from hot gas and X-ray binaries originating 
in a starburst), and model (B), a single-temperature (0.5-1 keV) thermal model + second 
power law. In addition, since emission lines are seen at 6.5-7.0 keV, narrow Gaussians at 6.4 
keV, 6.7 keV, and 7.0 keV, Ka emission lines from Fe at different ionized levels, were also 
included in the model fitting. Table 2 summarizes the best-fit models that were adopted. 
The Lx(2-10 keV) depends only weakly on the adopted models. We obtained EW6.4 ~ 340 
eV, N H ~ 3 x 10 23 cm" 2 , and Lx(2-10 keV) ~ 9 x 10 42 ergs s" 1 in both models. The Lx(2- 
10 keV) is about twice as large as the estimate based on ASCA data (Ueno et al. 1996). 
The absorption-corrected 2-10 keV X-ray to observed infrared luminosity ratio is Lx(2-10 
keV)/L IR ~ 5 x 10~ 3 . 

For PKS 1345+12, a second power law and a narrow Gaussian (6.4 keV Fe Ka line) 
were included, in addition to the main absorbed power-law component (model C). The X-ray 
spectrum is well fitted (Table 2). We obtained EW 6 . 4 ~ 130 eV, N H ~ 5 x 10 22 cm -2 , and 
Lx(2-10 keV) ~ 4 x 10 43 ergs s -1 . The Lx(2-10 keV) we obtain is about 30% higher than 
that derived from ASCA data (Imanishi & Ueno 1999b). Lx(2-10 keV)/L m is ~6 x 10~ 3 . 

For IRAS 05189—2524, model (D), a low temperature (0.5-1 keV) thermal model + 
an even lower temperature (<0.3 keV) second thermal component, and model (E), a single- 
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temperature (<1 keV) thermal model + second power-law component, were used, primarily 
to account for the soft X-ray emission. A narrow Gaussian for the Fe Ka line at 6.4 keV was 
also added, although this line is not as clear as in Mrk 463. Table 2 summarizes the adopted 
best-fit models. We obtained EW 6 . 4 ~ 30 eV, N H ~ 6 x 10 22 cm" 2 , and Lx(2-10 keV) 
~ 2 x 10 43 ergs s _1 . The Lx(2-10 keV) agrees with the ASCA and BeppoSAX estimates 
(Severgnini et al. 2001). Lx(2-10 keV)/L IR is ~4 x 10~ 3 . 

IRAS 07598+6508 shows no clear signature for an absorbed hard component. The 
spectrum was fitted with a primary power law and a narrow Gaussian for the Fe Ka 6.4 
keV line (Table 2). Although no strong Fe Ka line is evident at 6.4 keV, EWg.4 could not 
be constrained strongly (EW6.4 < 13 keV), due to poor photon statistics. Lx(2-10 keV) and 
Lx(2-10 keV)/L IR are estimated to be ~8 x 10 41 ergs s' 1 and ~9 x 10~ 5 , respectively. BAL 
quasars are known to be weaker X-ray emitters than non-BAL quasars, due to absorption 
(Green et al. 2001). The substantially smaller Lx(2-10 keV)/L m in IRAS 07598+6508 (BAL 
quasar) than the other three ULIRGs may be related to the BAL phenomena. The best fit 
value for the power law index T ~ 2.9 (Table 2) is larger than the canonical value for AGNs 
(r ~ 1.8-2.0). The large T value may be caused by the presence of additional soft X-ray 
emission which can steepen the apparent X-ray spectrum. However, the limited photon 
statistics hamper more detailed investigations of the spectrum. 



5. Discussion 

5.1. The 2—10 keV to Broad Line Luminosity Ratio 

For Mrk 463, PKS 1345+12, and IRAS 05189-2524, the modest EW 6 . 4 (< several x 
100 eV) and clear signatures for absorption with N H < several x 10 23 cm' 2 suggest that the 
detected hard component is direct emission from an AGN behind Compton-thin absorbing 
material (Awaki et al. 1991; Krolik, Madau, & Zycki 1994; Ghisellini, Haardt, & Matt 1994), 
and so the absorption-corrected 2-10 keV luminosities, Lx(2-10 keV), estimated in §4, can 
be taken as the intrinsic X-ray luminosity of the AGN. Therefore, the 2-10 keV to broad-line 
luminosity ratios for the AGNs in these ULIRGs can be compared with those for moderately 
infrared- luminous type-1 AGNs (Ward et al. 1988). For IRAS 07598+6508, however, the 
constraint on EW6.4 (<13 keV) is so loose that it is unclear whether the observed 2-10 keV 
flux is a direct component behind Compton-thin absorbing material (EW 6 . 4 < several x 100 
eV). The extremely small Lx(2-10 keV)/Li R ratio and small N H (~8 x 10 20 cm~ 2 ) suggest 
that the observed 2-10 keV emission is a scattered/reflected component behind Compton- 
thick absorbing material. This source will be excluded from the following discussions. 
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The correlation between 2-10 keV and broad emission lines established for type-1 AGNs 
(Ward et al. 1988) is for the broad component of an optical Ha (A rest = 0.65 //m) emission 
line. Thus, we need to convert the luminosity of a broad near-infrared emission line (Pact for 
PKS 1345+12 and IRAS 05189-2524 2 ; Veilleux et al. 1997b,1999b, and Pa/? for Mrk 463; 
Veilleux et al. 1997a) to that of broad Ha. For emission lines from the broad-line regions in 
an AGN, the case-B assumption is not applicable to low-level transition lines, such as Ha, 
due to collisional effects (Osterbrock 1989). For the broad Ha to Paa luminosity ratio, the 
value of 20 observed in dust-unobscured type-1 AGNs (Hill et al. 1996) was adopted. Since 
the deviation of this ratio from case-B (~10) is only a factor of ~2, the choice of this ratio 
will not affect our main conclusions significantly. For high-transition broad emission lines, 
such as Paa and Pa/3, case-B is applicable (Rhee & Larkin 2000). Thus, a broad Paa to 
Pa/3 luminosity ratio of 1.9 (case-B) was adopted. 

Absorption- corrected 2-10 keV and extinction- uncorrected broad emission-line luminosi- 
ties for the three ULIRGs are plotted in Figure 3 (solid squares), where the estimated 
broad Ha emission-line luminosities are ~5, ~3, and ~1 x 10 43 ergs s -1 for Mrk 463, PKS 
1345+12, and IRAS 05189-2524, respectively. The type-1 AGNs studied by Ward et al. 
(1988) are plotted for comparison (open circles). In addition to the three ULIRGs, the 
other two Seyfert-2 ULIRGs with detectable broad near-infrared Paa emission lines, IRAS 
20460+1925 and IRAS 23060+0505 (Veilleux et al. 1997b), have available ASCA 2-10 keV 
spectra (Ogasaka et al. 1997; Brandt et al. 1997). Since their observed 2-10 keV fluxes are 
sufficiently high (>10~ 12 ergs s _1 cm -2 ), strong constraints were obtained for EW 6 . 4 (260^^7 
eV for IRAS 20460+1925; Ogasaka et al. 1997, <290 eV for IRAS 23060+0505; Brandt et 
al. 1997). The detected 2-10 keV fluxes show signatures for absorption with Nh = 2-8 
x 10 22 cm -2 , so that the detected 2-10 keV flux is taken as a direct component from an 
AGN behind Compton-thin absorbing material. For IRAS 20460+1925, Lx(2-10 keV) is 
~1 x 10 44 ergs s- 1 (Ogasaka et al. 1997) and Lx(2-10 keV)/L IR is ~1 x 10~ 2 . For IRAS 
23060+0505, Lx(2-10 keV) is ~2 x 10 44 ergs s" 1 (Brandt et al. 1997) and Lx(2-10 keV)/L IR 
is ~2 x 10~ 2 . Broad Ha emission-line luminosities were derived from Paa data (Veilleux 
et al. 1997b) using the method described above and estimated to be ~2 and ~1 x 10 44 
ergs s" 1 for IRAS 20460+1925 and IRAS 23060+0505, respectively. These two ULIRGs are 
also plotted in Fig. 3 (solid squares). The five Seyfert-2 ULIRGs show absorption- corrected 
2-10 keV to extinction-uncorrected broad line luminosity ratios roughly an order of magni- 
tude smaller than those in dust-unobscured type-1 AGNs with modest infrared luminosities. 
While no significant dust extinction is expected for broad emission lines in type-1 AGNs, 



2 In IRAS 05189—2524, broad components were detected at Ha, Pa/3, and Paa (see § 2). We use the Paa 
line, because dust extinction is the lowest. 
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some degree of dust extinction may be present for the five Seyfert-2 ULIRGs. For these 
five ULIRGs, since near-infrared Pact or Pa/3 emission lines were used to estimate the broad 
Ha luminosities, the effects of dust extinction are expected to be reduced. If some amount 
of dust extinction is present for the broad near-infrared emission lines and its correction is 
applied, these five ULIRGs will move to the right in Fig. 3, making the deviation of these 
ULIRGs even larger. Even using new XMM and Chandra data, we confirmed systematically 
smaller 2-10 keV to broad line luminosity ratios in ULIRGs, originally suggested by Imanishi 
& Ueno (1999a). 

It is unlikely that the systematic difference is caused by the uncertainty of the luminosity 
conversion from broad Paa or Pa/3 to Ha for the following two reasons. First, in a number 
of moderately infrared-luminous type-2 AGNs, the 2-10 keV to broad-Ha luminosity ratios 
derived from observed broad near-infrared line luminosities agree with those for unobscured 
type-1 AGNs (Imanishi & Ueno 1999a). Second, for IRAS 05189-2524, the intrinsic broad 
Ha luminosity estimated from the observed broad Paa luminosity (Lh q ~ 1 x 10 43 ergs s _1 ; 
see previous paragraph) is in good agreement with that derived directly from the broad Ha 
luminosity (Young et al. 1996) after correction for dust extinction estimated by Severgnini 
et al. (2001) is applied (Lh« ~ 1 x 10 43 ergs s _1 ). Both of these results suggest that (1) 
the luminosity conversion from broad near-infrared emission lines to broad Ha is reasonable, 
and that (2) at least for some sources with detectable broad near-infrared emission lines, 
dust extinction is not significant in the near-infrared domain. Thus, the significant deviation 
of ULIRGs in Fig. 3 is believed to be real. 

5.2. X-ray Underluminous or Broad Line Over luminous? 

Estimating the fractional AGN contribution to the large infrared luminosities of ULIRGs 
is of particular importance to the study of their energy sources. It is not clear which of the 
two AGN indicators, 2-10 keV emission or broad optical/near-infrared emission lines, is a 
better tracer of this fraction. The 2-10 keV luminosity depends primarily on the number 
of hot electrons in the close vicinity of an accretion disk around a central supermassive 
blackhole in an AGN and the number of incident UV photons to be Compton up-scattered. 
The broad emission-line luminosity is roughly proportional to the product of the number of 
UV photons from the central AGN and the covering factor of broad-line gas clouds. The 
infrared luminosity of a ULIRG powered by an AGN is determined mainly by the number 
of UV photons from the AGN and by the covering factor of dust. In the nuclear region 
of a ULIRG, a large amount of gas and dust is available (Sanders & Mirabel 1996), which 
may increase the covering factors both of dust and broad-line gas clouds around the AGN. 
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If this is the case, broad emission lines are a better indicator of the AGN's contribution to 
the infrared luminosity, and AGNs in ULIRGs are underluminous in 2-10 keV X-rays. 

In fact, in a number of ULIRGs, although an AGN is consistently suggested to be 
the primary energy source by some independent methods, only energy diagnostics based 
on the 2-10 keV luminosity suggest that the AGN is energetically insufficient. Examples 
include IRAS 05189-2524 (Veilleux et al. 1999b; Soifer et al. 2000; Imanishi & Dudley 2000; 
Severgnini et al. 2001) and Mrk 463 (Ueno et al. 1996; Genzel et al. 1998; Imanishi 2002). 
Thus, it is more likely that the AGNs in ULIRGs are underluminous in X-rays with respect 
to their overall spectral energy distributions, than that their broad lines are overluminous. 

If X-ray underluminosity (relative to UV) is an intrinsic property of ULIRGs, one possi- 
ble explanation is a high mass accretion rate onto a central supermassive black hole (Bechtold 
et al. 2003). This is plausible in ULIRGs, given the high nuclear concentration of gas and 
dust (Sanders & Mirabel 1996). 

Alternatively, the apparent X-ray underluminosity of AGNs in ULIRGs may be caused 
by our misidentification of the origin of the detected 2-10 keV emission from ULIRGs. 
Using reasonable assumptions for an AGN, the scattered/reflected X-ray component of AGN 
emission behind Compton-thick absorbing material is expected to show a very large EW6.4 
with >1 keV and only weak absorption (Matt, Guainazzi, & Maiolino 2003). A modest EW 6 . 4 
and clear Compton-thin absorption signature, as were found in the five ULIRGs, make it 
most plausible that we are observing direct X-ray emission from an AGN behind Compton- 
thin absorbing material (Awaki et al. 1991; Krolik, Madau, & Zycki 1994; Ghisellini, Haardt, 
& Matt 1994). However, Thomson scattering by highly ionized gas produces a small EW6.4. 
If the scattered component suffers some degree of absorption, both the modest EW6.4 and 
clear signature for Compton-thin absorption can be explained (Maloney & Reynolds 2000). 
In this scenario, the intrinsic 2-10 keV luminosity of an AGN could be substantially larger, 
depending on the scattering efficiency Since such X-ray spectra are usually interpreted as 
Compton-thin, and the absorption-corrected 2-10 keV luminosity is determined based on 
this assumption, we will underestimate the true Lx(2-10 keV). 

A third possibility which could explain the observed X-ray underluminosity is partial 
Compton-thick X-ray absorption. If the X-ray-absorbing material in front of an AGN consists 
of small clumps of Compton-thick and thin clouds, whose size scales are substantially smaller 
than the size of the X-ray-emitting region of the AGN, then the 2-10 keV spectrum will 
suggest Compton-thin absorption. The actual Lx(2-10 keV) could be significantly larger 
than a simple estimate assuming uniform Compton-thin X-ray absorption. The covering 
factor of Compton-thick clouds in front of an AGN along our line-of-sight determines the 
degree of our underestimate for Lx(2-10 keV), but an X-ray spectrum at >10 keV is required 
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to determine this (Braito et al. 2003b). 

Whether the X-ray underluminosity of AGNs in ULIRGs is intrinsic or caused by our 
underestimate of Lx(2-10 keV), the use of 2-10 keV observations to estimate AGNs' con- 
tribution to the infrared luminosities of ULIRGs will be quite misleading, unless the X-ray 
underluminosity is properly taken into account. 

The discovery of underluminous X-ray emission relative to broad emission lines from 
AGNs is only for ULIRGs with modest dust obscuration along our sightlines, because in order 
for broad optical/near- infrared emission lines to be detected, dust obscuration toward AGNs 
must be relatively low. The putative AGNs in the majority of ULIRGs are so highly dust- 
obscured that broad emission lines are not detectable (Veilleux et al. 1999b), and so we have 
no information on the broad line luminosities. It is unclear whether X-ray underluminosity 
is a general property of AGNs in the majority of ULIRGs, and there is currently no way to 
determine this. However, when using 2-10 keV data to determine the energy contribution of 
AGNs in ULIRGs, it is important to consider that the AGN contribution may be significantly 
underestimated. 



6. Summary 

We estimated the absorption-corrected 2-10 keV X-ray luminosities of AGNs in mod- 
erately dust-obscured ULIRGs, and confirmed a previous finding that absorption- corrected 
2-10 keV to extinction- uncorrected broad emission-line luminosity ratios in ULIRGs are 
roughly an order of magnitude smaller than those of type-1 AGNs with modest infrared 
luminosities. We have shown that the smaller ratios are more likely to be due to 2-10 keV 
X-ray underluminosity with respect to the overall spectral energy distributions, rather than 
higher broad-emission-line luminosity. The lower X-ray luminosities of AGNs in ULIRGs 
should be taken into account in estimating AGNs' contribution to the energy output from 
ULIRGs based on X-ray data. 
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of the SIMBAD database, operated at CDS, Strasbourg, France, and of the NASA/IPAC 
Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, Cal- 
ifornia Institute of Technology, under contract with the National Aeronautics and Space 
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Table 1. Details of Observed ULIRGs. 



Object 


z 


Lir 


Satellite 


Date 




Net Exposure 






[ergs s -1 ] 




[UT] 




[ksec] 


(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


Mrk 463 


0.051 


45.3 1 


XMM 


2001 Dec 


22 


21 (PN), 25 (MOS) 


PKS 1345+12 


0.122 


45.8 


Chandra 


2000 Feb 


24 


19 


IRAS 05189-2524 


0.042 


45.7 


XMM 


2001 Mar 


17 


7 (PN), 10 (MOS) 


IRAS 07598+6508 


0.149 


46.0 


XMM 


2001 Oct 


25 


13 (PN), 19 (MOS) 



Note. Column (1): Object. Column (2): Redshift. Column (3): Logarithm of 
infrared (8—1000 /im) luminosity in ergs s" _1 calculated with L m = 2.1 x 10 39 x D(Mpc) 2 
x (13.48 x /12 + 5.16 x / 25 + 2.58 x / 60 + /ioo) ergs s" 1 (Sanders & Mirabel 1996), where 
fi2 , f25j f60j and boo ar e IRAS FSC fluxes at 12/xm, 25/um, 60/um, and lOO^um, respectively. 
Column (4): X-ray satellite used for the observation. Column (5): Observing date in UT. 
Column (6): Effective exposure time in ksec after discarding high background intervals. 

x Mrk 463 is not a ULIRG in the strict sense, but is included, since its infrared luminosity 
is close to that of a ULIRG. 
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Table 2. Results of X-ray Spectral Fitting. 



Object 


Model 1 


X 2 /d.o.f. 


N H 


L x (2-10keV) 


EW 6 . 4 








[10 22 cm" 2 ] 


[10 42 ergs s _1 ] 


[eV] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Mrk 463 


A 


231.6/202 


33lg 


9.2 


3401™ 




B 


230.4/203 


32+3 


9.1 


3401™ 


PKS 1345+12 


C 


49.2/62 


4 ^+0-5 
4 -°-0.4 


38 


130+130 


IRAS 05189-2524 


D 


225.9/232 


5.7+0.3 


17 


30H|8 




E 


233.2/233 


5.8+0.4 


16 


301™ 


IRAS 07598+6508 


F 


25.1/23 


o-o8±8:S8 


0.8 


< 13000 



Note. — Column (1): Object. Column (2): Adopted best-fit models. Column (3): Reduced x 2 
values. Column (4): Absorption for the main power law component. The uncertainty is at 90% 
confidence level for one parameter of interest (A% 2 = 2.7) throughout this Table. Column (5): 
Absorption-corrected 2-10 keV X-ray luminosity for the main power law component, estimated 
based on the assumption of uniform Compton-thin absorption. Column (6): Equivalent width of 
the 6.4 keV Fe Ka emission line relative to the main power law component. 

1 Model A : absorbed power law (r = 1.8 fixed) + thermal (kT = 0.6+0.1 keV) + 
thermal (kT = 10 keV fixed) + narrow Gaussian (6.4 keV Fe Ka) + narrow Gaussian 
(6.7 keV, EW = 30±|g eV) + narrow Gaussian (7.0 keV, EW = 120±Jg° eV) 

Model B : absorbed power law (r = 1.8 fixed) + thermal (kT = 0.6+0.1 keV) 
+ second power law (r = 1.8 fixed; ~5% of the main power law component in 
normalization) + narrow Gaussian (6.4 keV Fe Ka) + narrow Gaussian (6.7 keV, EW 
= 40l|g eV) + narrow Gaussian (7.0 keV, EW = 130+80 eV) 

Model C : absorbed power law (r = 1.8 fixed) + second power law (r = 1.8 
fixed; ~3% of the main power law component in normalization) + narrow Gaussian 
(6.4 keV Fe Ka) 

Model D : absorbed power law (r = 1.8 fixed) + thermal (kT = 0.7^ keV) 
+ thermal (kT = 0.11+0.02 keV) + narrow Gaussian (6.4 keV Fe Ka) 

Model E : absorbed power law (r = 1.8 fixed) + thermal (kT = 0.3+0.1 keV) 
+ second power law (r = 1.8 fixed; ~1% of the main power law component in 
normalization) + narrow Gaussian (6.4 keV Fe Ka) 

Model F : absorbed power law (r = 2.91q'|) + narrow Gaussian (6.4 keV Fe 
Ka) 



Fig. 1. — A Chandra 0.5-8 keV X-ray image of PKS 1345+12, with infrared 2.2-^xm contours 
obtained from archival HST NICMOS data overplotted. North is up and east is to the left. 
The field of view is 12 x 12 arcsec 2 . The contours from the lowest one are 12.5%, 25%, 
50%, and 100% of the highest level at the western nucleus. The coordinate of the X-ray 
emission peak in this Chandra image is J2000.0 [13:47:33.37, 12:17:24.3], which is slightly 
displaced from the 2.2-/zm peak of the western nucleus estimated from the NICMOS archival 
data (J2000.0[13:47:33.33, 12:17:24.0]). Evans et al. (1999) made more accurate astrometry 
of this galaxy, and estimated the coordinate of the 2.2 /mi peak of the western nucleus 
to be J2000.0[13:47:33.38, 12:17:24.4]. Considering the absolute positional uncertainty of 
Chandra (0.6 arcsec for 90% error radius), we conclude that the X-ray emission comes from 
the western nucleus. 
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Fig. 2. — (a): An XMM spectrum of Mrk 463. Model A is overplotted as solid lines. In the 
upper panel, the higher plots are a PN spectrum (which contains the most counts), and the 
lower plots are MOS spectra. The lower panel shows residuals of data from the model, (b): 
A Chandra spectrum of PKS 1345+12. Model C is overplotted. (c): An XMM spectrum of 
IRAS 05189-2524. Model D is overplotted. (d): An XMM spectrum of IRAS 07598+6508. 
Model F is overplotted. 
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Fig. 3. — Absorption- corrected 2-10 keV X-ray luminosity (ordinate) and extinction- 
uncorrected broad Ha emission line luminosity (abscissa). Open circles are dust-unobscured 
type-1 AGNs with moderate infrared luminosity taken from Ward et al. (1988). Filled 
squares are the five Seyfert-2 ULIRGs studied in this paper, Mrk 463, PKS 1345+12, IRAS 
05189-2524, IRAS 20460+1925, and IRAS 23060+0505 (see text). The solid line is for log 
Lx(2-10 keV) = log L(broad Ha) + 1.255 (Imanishi & Ueno 1999a). In the dashed lines, 
Lx(2-10 keV)-to-L (broad Ha) ratio is a factor of 6 larger and smaller than the solid line. The 
only X-ray underluminous AGN (IRAS 05218-1212) originally found by Ward et al. (1988) 
is infrared selected by IRAS, if not a ULIRG. Since no information on EW 6 . 4 is available, 
the weak 2-10 keV X-ray emission might be caused by Compton-thick absorption. 



